UNCLASSIFIED 


*‘»428117 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  AlEXANDRIA.  VIRGINIA 


UNCLASSIFIED 

AOoSOillOSI 


■onci:  Mbea  tanramtat  or  othar  dra«laci>  apael- 
fleatlGoa  or  othar  data  are  uaad  for  aay  puxpoaa 
other  thaa  In  conneetion  vlth  a  deflaltaljr  related 
aoraiuaent  proeuranent  operation,  the  U.  8. 
flmenwant  thereby  laeurr  no  ntapooslbillty,  aor  any 
ObUcatlon  ehataoerar}  and  the  fleet  that  the  Oovam- 
■ent  nay  have  foaailated,  fUnii^ied,  or  la  asy  wy 
ct^plied  the  aald  drawlacB,  epeelfleationa,  or  other 
data  le  not  to  be  repurdad  by  laplleattsa  or  other- 
vlae  ae  la  aay  nanaer  Ueeaelac  the  holder  or  any 
other  pereoo  or  corporation,  or  eoareyiag  aay  ridhte 
or  pemiaaioa  to  aaaufiaetvre,  uae  or  aall  any 
pateated  iarention  that  aay  ia  any  nay  be  related 
thereto. 


SPFCIAL  ASPECTS  OF  ENVIRONHENT  RESULTING 
FROM  VARIOUS  KINDS  OF  NUCLEAR  WARS 

PART  1 1  Annex  1 1 1 

Application  of  Input-Outpgt  Analysis 
to  a  Homeostatic  Ecosvitem 

Robert  U.  Ayres 

HI-303-RH/A  III  January  8,  1964 


Prepared  under  Contract  No.  0C0*0S-62-2l8 
Department  of  Defense,  Office  of  Civil  Defense 
OCD  Task  Number  3$11A 


This  report  has  been  reviewed  In  the  Office  of  Civil 
Defense  and  approved  for  publication.  Approval  does 
not  signify  that  the  contents  necessarily  reflect  the 
views  and  policies  of  the  Office  of  Civil  Defense. 


Qualified  requesters  may  obtain 
copies  of  this  report  from  DOC. 


HUDSON  INSTITUTE,  INC. 
Quaker  Ridge  Road 
Ha  rmon-on- Hud son 
New  York 


PUBLtCATIONS  IN  THIS  SERIES 


SPECIAL  ASPECTS  OF  ENVIRONHENT  RESULTING  FROM  VARIOUS  KINDS  OF  NUCLEAR  WARS 

i 


Part  I  Including  Annexes  l-V  and  Appendix  I 

/  HI-2l*3-RR  Junes,  1963  ..R.U.  Ayres 

Appendices  2  and  3:  HISTORICAL  EXAMPLES  OF  ECOLOGICAL  DISASTER 

HI-243-RR/A  2-3  June  5,  1963.,... . J.  Ingersoll 


Part  II  including  Annexes  1,11:  THE  USE  OF  SCENARIOS  FOR  EVALUATING 

POSTATTACK  DISUTILITIES 

HI-303-RR  January  8,  1964  . R.U.  A, res 

Annex  111:  APPLICATION  OF  INPUT-OUTPUT  ANALYSIS  TO  A  "OHEOSTATIC 

ECOSYSTEM 

HI-303-RR/A  III  January  8, .  1964.. . . . . .  ,R.U.  Ayres 

Annex  IV:  EFFECTS  OF  A  THERMONUCLEAR  WAR  ON  THE  WEATHER  AND  CLIMATE 

HI-3O3-RR/A  IV  January  8,  1964. . R.U.  Ayres 

Appendixes  1,2:  HISTORICAL  EXAMPLES  OF  ECOLOGICAL  DISASTER  (II) 

HI-303-RP'A  1-2  January  8,  1964 . 


J,  Ingersoll 


HI-303-RR  I 

ANNEX  til 

APPLICATION  OF  INPUT-OUrPUT  ANALYSIS 
TO  A  HOMEOSTATIC  ECOSYSTEM 


Input-output  analysis  1$  a  technique  which  has  been  rather  highly 
developed  in  economics,*  It  has  been  quite  successful  In  making  certain 
sorts  of  economic  predictions,  though  it  is  not  and  does  not  purport  to 
be  a  complete  treatment  of  that  subject. 

In  ecology  it  seems  reasonable  to  conjecture  that  input-outpu, 
analysis  might  also  be  appropriate  in  treating  rertain  type':  of  ques¬ 
tions.  In  particular,  one  would  like  to  ask  such  things  as,  "What 
level  of  insect  feeding  on  crops  would  be  iikely  to  result  if  the  insec¬ 
tivorous  birds  were  driven  away,  all  other  conditions  remaining  unchanged?" 
or  a  more  general  version,  "What  is  the  likely  consequence  of  such-and- 
such  a  selective  pesticide  removing  some  species  and  not  others?" 

Such  questions  not  only  have  many  obvious  peacetime  agricultural 
appi icat 'ons ,  but  also  bear  on  the  design  of  sophisticated  life-support 
systems  and  even  one  of  the  most  difficult  of  all  the  problems  involved 
in  a  study  of  a  hypothetical  postattack  environment,  namely  the  effect 
of  a  selective  toxic  agent  (e.g.  radioactive  fallout)  on  the  "balance 
of  nature." 

The  possibility  of  using  this  type  of  analysis  in  the  ecosphere 

I  I 

depends  on  the  fac  :  that  biological  production  is  in  fact  limited  by 
the  availability  o'  certain  essential  elements  in  thr  s',stsm,  i.e., 
energy,  water,  nitrogen,  phosphorus,  carbon,  etc.  Where  other  factors 
are  dominant,  input-output  analysis  would  cleariy  be  inapplicable.  For 
example,  the  rate  of  growth  of  bacteria  in  a  culture  or  the  nature  of 
the  biological  succession  in  an  "old  field"  are  determined  by  dif- 
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ferent  cons tderat ions.  These  are  nonequi I tbrlum  problems. 

Generally  speaking.  Input-output  would  seem  to  be  most  appropriate 
where  a  homeostatic  "steady-state"  has  been  reached,  either  naturally 
(e.g.  climax,  vegatatlon)  or  as  a  result  of  the  institutionalized 
permanent  Interference  of  man.  By  steady-state  in  this  context  we  mean 
that  a  pattern  is  repeated  year  after  year,  on  the  average.  The  annual 
cycle  of  growth  can  perhaps  be  taken  Into  account  adequately  by  assuming 
a  simple  sinusoidal  oscillation  of  the  relevant  variables  with  a  one- 
year  period.  Blennl#!,  triennial  or  quaJrjennial  cycles  (e.g.  crop 
rotation  schemes)  could  perhaps  be  handled  in  a  similar  way,  but  a 
simpler  and  perhaps  equally  valid  procedure  would  be  to  average  over  a 
relatively  homogeneous  area  (perhaps  a  county)  and  assume  the  distri¬ 
bution-percentage  of  land  devoted  to  different  crops— Is  constant  from 
year  to  year  even  though  the  actual  crops  on  a  particular  land-parcel 
may  change. 

In  comparison,  most  other  variations  with  time  can  either  be  averaged 
out  or  neglected,  at  least  for  an  agricuTttirai  community  or  a  complex 
temperate-zone  ecosystem  such  as  a  mixed  forest.  This  would  apply  to 
the  day-night  ''ariation— a  kind  of  high  frequency  "hufr."-*and  to  random 
short-term  deviations  from  the  average  die  to  weather  fluctuations— 
analogous  to  "noise".  Long-term  climatic  fluctuations  and  low-fre¬ 
quency  cyclic  variations  (e.g.  the  sunspot  cycle  and  others)  are  some¬ 
what  harder  to  take  into  account.  We  will  assume  for  the  present  that 
these  variations  have  small  amplitudes  and  can  be  neglected,  although 
in  some  well  known  instances  this  assumption  would  clearly  be  faulty. 

An  alternative  procedure  will  be  discussed  which  requires  no  assump¬ 
tions  regarding  the  time-behavior  of  the  model  ecosystem  and  can  be 
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applied.  In  principle,  even  to  non-homeostatlc  situations,  at  the  cost 
of  considerable  labor. 

$2  Development  of  the  Model 

Bases  of  Transactions:  We  focus  attention  on  one  of  the  elements 
being  exchanged  between  various  components  of  a  biological  system,  A 
number  of  such  elements  exist,  which  may  be  labeled 

• 

For  purposes  of  convenient  reference,  suppose  ^  represents  nitrogen; 

phosphorus;  sulfur,  etc.  Carbon,  water,  energy,  calcium  et  al . 
are  alt'  possible  choices.  In  any  particular  case,  our  choice  or  choices 
of  which  elements  are  critical  must  be  guided  by  prior  knowledge  of  the 
system. 

Components;  The  elementary  components  of  any  ecosystem  are  popu¬ 
lations  of  individual  species.  In  practice,  however,  much  broader  and 
more  Inclusive  groupings  are  necessary  if  the  problem  is  to  be  trac¬ 
table  (even  with  the  help  of  electronic  computers). 

The  optimal  division  of  the  system  into  lumped  classes  is  not 
altogether  well  defined.  Two  criteria  are  i.nportant,  neither  of  which 
necessarily  would  lead  to  the  most  aesthetically  "natural"  list  of  cate- 
gorles--which  tpnjuld  differ  f.ir  different  biologists,  but  would  probably 
be  based  in  some  sense  on  current  understanding  of  phylogeny.  It  Is 
vital  to  realize  that  for  an  Input-output  study  phyTogeny  is 
tially  irrelevant.  The  two  criteria  are: 

(!)  If  the  model  Is  to  answer  questions  about  some  particular 
class,  that  class  cannot  be  lunped  Into  a  larger  group. 

(2)  The  categories  must  be  chosen  In  such  a  way  that  it  makes 
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sencc  to  describe  an  Interaction  between  them  in  terms  of  exchange  of 
one  of  the  critical  substances,  e.g. 

Point  (I)  is  best  illustrated  by  an  example:  if  we  are  especially 
interested  in  the  importance  of  bats  in  controlling  night-flying  insects, 
then  both  must  appear  explicitly,  i.e.,  we  cannot  lump  bats  with  mam¬ 
mals  or  vertebrates,  and  night-flying  insects  must  explicitly  be  distin¬ 
guished  from  other  types. 

Point  (2)  is  less  easy  to  explain  precisely.  The  best  illustration 
may  be  to  give  an  example  of  an  interaction  which  is  not  expressible  in 
terms  of  exchanging  any  simple  substance  which  has  been  mentioned  so 
far:  growth  of  the  understory  in  a  forest  is  determined  primarily  by 
the  amount  of  light  reaching  it  through  the  canopy.  A  division  into 
"understory  species"  and 'tanopy  species"  would  not  satisfy  the  criterion 
unless  represented  light  —  in  which  case  most  other  interactions  could 
not  be  expressed  satisfactorily. 

Assuming  an  appropriate  division  for  the  model  ecosystem  can  be 
achieved,  each  class  of  components  contains-  a  quantity  of  the  basic 
exchangeable  substances  ^ #  5  •  •  •  there  are  M  lumped  classes 
In  all,  these  quantities  at  any  time  are  given  by  the  variables 


For  simplicity  we  shall  discard  the  ^  subscripts 

Icwlng.  Since  5.7-5  circulate  constantly  throughout 
between  any  two  arbitrarily  chosen  classes  of  components 


in  the  fol- 
the  system, 
there  Is 


likely  to  be  a  net  flow  of  the  exchangeable  substance, in  one  direction 

of  the  other.  Let  us  consider  the  flow  of  ^  between  the  and 

>'^e  use  the  v/ord  "contains"  in  lieu  of  a  more  accurate  but  cuMiberso"ie 
terminolog”.  The  operational  definitions  of  the  A|^'s  are  discussed  in  ^4. 


classes  of  components.  This  transfer  is  characterized  by  a  coefficient 


or  "matrix-element"  which  describes  a  channel  having  a  uni-direc¬ 

tional  character,  allowing  ^  to  move  •.'•ly  from  J  to  k.  A  channel  in 
the  reverse  direction  Is  described  by  and,  in  general,  *ijk. 

Here  rfjk  represents  the  fraction  of  all  ^  *°"tgo  from  the  lumped- 

class  going  to  the  k^^  lumped-class  in  a  time  increment  St.  If  we 
allow  the  possibility  of  internal  transactions  within  a  lumped-class, 
then  the  following  relations  must  hold  true  if  all  interactions  are 
taken  into  account: 


where  4  0  in  general. 

Open  and  Closed  Systems;  Reservoirs:  Equation  (I)  characterizes 
a  closed  system  in  which  the  exchangeable  substance  ^  is  conserved. 
Most  ecosystems  are  (conceptually)  open  systems,  however.  That  is, 
there  is  a  constant  shuttle  back  and  forth  between  the  (open)  system 
and  an  external  reservoir  of  The  reservoir  may  be  thought  of,  in 
general,  as  the  soil,  rocks,  and  the  atmosphere,  i.e.,  the  inorganic 
world.  It  may  be  convenient  to  include  the  works  of  man,  as  well, 
unless  production-consumption  relationships  are  being  investigated. 

In  an  open  system  equation  (1)  will  not  hold,  instead,  there  will 
be  a  system  n*  inequalities 


k-i 


I  for  all  j 


lb 


The  N'«N-I  notation  indicates  that  one  component  (the  reservoir)  has 
been  removed  from  consideration,  whence  transactions  to  and  from  the 
resenoi  r  are  omitted.  Henceforth  N  will  be  used  indiscriminately  and 
distinctions  (if  any)  between  closed/open  systems  will  be  explicit. 


A 
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The  matrix  A  describes  the  transfer  of  in  the  system  in 
terms  of  the  "sources"  and  "sinks"  of  ^  .  One  immediately  obtains  a 


set  of  equations: 


<<ii,(A’  -  )  for  al  I  k 


T—  -  ••  i  k  \ry'  -  n  /  i  vi  a  i  f  t . 

t  j=l  J 

which  expresses  the  fact  that,  in  an  increment  of  time  ^t  any  excess 
or  deficit  in  the  quantity  of  ^  contained  in  the  k^*'  lumped-class  of 
components,  ^A|^(^)  must  have  originated  from  the  other' N- 1  components, 
or  reappeared  in  them.  Note  that  internal  transactions  do  not 

contribute  either  v/ay. 


Continuous  Cyclic  Model:  In  the  limit  as  gt  >0,  the  equation 


becomes 


\  ^  <ljk(Aj  -  Aj^)  -  (  ^  3. 

j-l  j»1 

where  A|^  is  the  time  derivative  of  A|^  and  the  's  are  now  Instantaneous 
rates.  The  most  generai  time  behavior  of  the  system  would  be  given  by 
a  Fourier  integral:  ^ 

*k  *k  e'*^*  AM 

-oo 

One  thing  we  know  about  the  system  from  the  outset  is  that  there 
exists  a  single  overwhelmingly  don;inant  frequency  CJq  with  a  period  of 
one  year,  which  expresses  the  most  salient  fact  about  the  time  behavior 
of  the  ecosystem:  namely,  the  annual  cycle.  The  best  approximate  tjme- 
function  to  pick  depends  on  the  climate  and  on  the  life  cycles  of  the 

organisms  making  up  each  component.  In  the  tropics  the  sinusoidal . 

variation  of  a  typical  A|^  is  very  slight,  being  superimposed  on  a  non¬ 
zero  average  value,  e.g.  fig.  (i): 
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4b. 


The  time-function  corresponding  to  fig.  (I)  might  be  of  the  form 
A(j(t)"  +  cosdJo^^  ^ 

In  a  temperate  climate,  the  curves  for  annuals  and  perennials  (both 
plant  and  animal)  are  quite  different.  The  curve  for  most  perennial', 
resembles  fig.  (I),  k-dilte  that  for  annual  species  Is  more  like  fig.  (it); 
e.g.  It  is  essentially  zero  during  the  winter  months  and  rises  sharply 
during  the  growing  season,  which  spans  the  period  between  frosts.  In 
the  mid-cont irental  U.S,  at  roughly  ktfi  N.  Latitude  the  growing  season 
is  quite  close  to  l80  days  or  half  of  tko  year.  For  an  ecosystem  con¬ 
sisting  largely  of  annual  species  at  such  a  latitude,  the  teme-dependence 
should  be  given  by  a  function  such  as 

A(j(t)«  A(^|C|5  +  cos  6^<,(t  +  A  k)0kt*j]'*’  — 
where  C^ls  a  constant  depending  on  climate,  (C  ^  0  at  k0°  N.  Latitude 
in  the  midvyestsrn  U.S.)y  ^  Is  a  very  small  positive  constant,  and 
^(0  Is  zero  when  C^^  +  cos  ^^^(t  +  \|j)  ^  0,  and  unity  otherwise. 

*  For  either  (4b)  or  (4c)  It  Is  trivial  to  show  that  the  integral  of 
A  over  a  full  c^le  vanishes,  viz. 

-  ^  Ak(t)  dt  -  Ak(T)  -  Ak(0)  -  0  5 

where  T  •  In  an  agricultural  community  in  the  temperate  zone 

(e.g.  Illinois)  the  choice  (4c)  with  C|(  •  0,  seems  most  applicable.  We 
shall  assune  this  In  the  following  for  simplicity. 
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Stmt !«rly ,  we  take 


.  0  ■ «( j )[« 


+  cos  <k»  (t 


where  and  Tj)^  unknown  phase  constants. 

Note  that  ft  Q  ■©  where  ©  Is  unity  over  the  range  where  both 
k  j  k  s  ''s 

e,  and  0j|^  are  unity,  and  zero  elsewhere.  Integrating  over  the  cycle, 

-  ■  I  ['  -  ‘I 


where  is  equal  to  unity  when  the  term  in  brackets  i  -  positive, 

and  zero  when  it  becomes  negative.  We  therefore  have,  for  the  case 
w  Cji^  a  o,  and  £  sufficiently  small  so  that  6^  can  be  neglected: 

°  ^jk  kAk  <=°*"o  (\  -Tjk)  -  \ 


8. 


J 


Note  that  Aj|^  always  has  the  same  Indices  as  the  argument  of  the 
accompanying  cosine  term,  so  the  subscripts  will  be  dropped  henceforth. 
This  is  a  set  of  linear  homogeneous  equations,  if  we  treat  A|...A^  as 
unknowns.  Nontrivial  independent  solutions  exist  only  If  the  deter¬ 
minant  of  the  matrix  vanishes, 


Oat  C<  «  0  where 


cos 


<rjk-\)-  s 


jk 


9. 


If  one  of  the  A's  is  fixed  arbitrarily  the  remaining  ones  can  be 
determined  by  straightforward  methods.  Altogether  there  are  N-I 
Independent  relations  between  the  A's  andCC’s  and  one  involving  only 
the  fit's  (equ.  9). 

It  should  be  noted  that  the  continuous  sinusoidal  model,  as  we 
have  derived  it,  cannot  be  applied  without  modification  to  a  closed 


sy&tTn,  though  it  Is  consistent  with  the  assumption  of  an  open  system. 
This  difficulty  arises  because  equation  (1)  and  equation  (^)— derived 
from  the  requirement  of  ^-conservation  in  a  closed  system— cannot  be 
satisfied  if  all  A's  and  A's  have  the  simple  sinusoidal  behavior 
assLined  in  (4)  and  (6).  A  more  complex  time  behavior  would  have  to  be 
assumed  for  at  least  one  of  the  A's  and  one  of  the  tljit's  (for  each  J), 
Discrete  Model :  A  different  model  is  derivab’e  from  the  basic 
equation  (2)  if  the  behavior  of  the  ecosystem  In  time  is  considered 
to  be  a  sequence  of  step-funct ions  Instead  of  a  smooth  variation. 
Equation  (2)  can  be  rewritten 

Mt  >'<yt)  -  \h)  \  ^  -  A^^b 

it  71  j-i 

Whence  (keeping  ^t  finite), 

A^(t  +  t)  -  Aj(t) 

or 

Ak'  -  tft  ^  djk  Aj  -  ^j^jkAj  10. 

j*.!  j-l'^ 

where  the  /^jk  represent  fractional  transfer  coefficients,  rather  than 
rates.  Note  that  (10)  can  be  iterated— thus  the  bahavior  of  the  system 
after  any  number  of  time  increments  is  given  by 

a(")  ./?<"-»•  •  .  a(°^  II. 

where  the  a(J)  are  vectors  and  the  are  matrices  of  rank  N  and 

■  A|^(t  +  n^t) 

Equations  (II)  supply  nN  equations  for  the  (n  +  l)N  different  Aj^'''^(s), 
which  is  sufficient  to  determine  all  Aj  (s)  for  m  ^  0,  assuming  the 
set  Aj are  given.  Not  that  if  the  ecosystem  undergoes  a  cyclic 
(annual  or  other)  variation,  the  Aj^  must  eventually  return  to  the  orl- 


12 


giriTl  set  of  values,  after  a  time  T  made  up  of  n  inert  •ents: 

-  A^(t  +  T)  -  A^^®^t) 
whence  ^ 

IT  -  1  >3. 

J-0  ^ 

(cycle) 

or,  in  words,  the  products  of  n  successive  /3-matrices  over  a  single 
period  In  a  cyclic  variation  must  be  equivalent  to  the  unit  matrix. 


Equation  (13)  is  the  analog,  in  the  discrete  model,  of  the  equation 
(7).  It  Is  equivalent  to 


^(n-l)  .(n-2) 

/S  ^ 

jl(i)  '>(i)'(2) 


(-.-!)“ 


r  -0Ji<k 


Equation  (14)  supplies  independent  conditions  between  the  nN^  dis¬ 
tinct  matrix  elements  y5jk^"^»  while  (II)  supplies,  in  this  Instance, 
nN  equations  relating  only  nN  different  Aj  values,  since  Aj 
Aj(®)  for  allj  by  the  cyclic  assumption.  Thus  a  certain  degree  of 
simplification  is  introduced  by  the  cyclic  condition,  but  r>  assunption 
of  a  particular  (e.g.  sinusoidal)  time  behavior  is  made  in  this  model. 
Moreover,  if  the  cyclic  condition  is  discarded,  the  discrete  model 
can  be  applied  to  Irreversible  or  non-cycllc  types  of  tine  behavior.* 


^3  Interpretation  and  Use  of  the  Models 

We  have  a  number  of  entities  (namely  l^'s,  A's  and  phase 
constants),  most  of  which  can,  in  principle,  be  measured  directly 


*Always  provided  the  ecosystem's  behavior  is  dominated  by  the  inter¬ 
change  of  some  substance  . 
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or  deduced  from  other  date.  Moreover  there  are  a  nuriber  of  additional 
equations  avcilabte  relating  the  A's  and  CL's  which  can  be  used  to 
reduce  the  nunber  of  independent  pieces  of  data  required. 

One  objective  ofour  analysis  is  to  examine  how  a  perturbation 
affects  the  ecosystem.  Many  types  of  perturbation  could  be  considered, 
but  the  simplest  one  which,  is  not  wildly  unreal  istic‘»  is  to  assume  that 
selective  damage  occurs  which  effectively  eliminates  one  or  more  of  the 
lunped  groups  of  components  without  directly  injuring  the  remainder. 

Our  assumption  about  the  simple  time-dependence  of  the  system  would  be 
invalid  until  the  system  had  "settled  down"  in  a  new  homeostatic  equili- 
briun.  We  can  ask,  however,  what  the  characteristics  of  the  new  system 
might  be,  i.e.,  the  values  of  the  A's  and  CL's  after  the  perturbation. 

Example;  Consider,  first,  the  continuous  cyclic  model  and  as  an 
il lustration  take  an  ecosystem  consisting  of  the  following  lumped- 
classes: 

Class  I.  Vegetation  (e.g. crops,  weeds) 

Decay  organisms 

2.  Plant-feeding  arthropods,  nematodes,  etc. 

3.  Predaceous  and  parasitic  arthropods,  R''  latodes,  etc. 

4.  Vertebrates  (birds,  mammals,  etc.) 

The  sample  system  Is  considered  open,  e.g.  the  total  N-content  is 
not  fixed,  but  shuttles  back  and  forth  between  the  system  and  a  nitrogen 

'^hls  type  of  perturbation  is  not  inconsistent  with  what  might  hap¬ 
pen  if  moderately  heavy  radioactive  fallout  were  to  cover  the  area  as  a 
result  of  a  nuclear  attack.  It  right  be  simulated  in  several  ways',  either 
by  very  active  hunting  or  trapping  of  birds  and  mannals  with  perhaps 
the  help  of  barriers  and  nets,  or  possibly  by  the  use  of  certain  poisons. 
It  is  known  that  birds  and  mammals  are,  on  the  average,  considerably 
more  radiosensitive  than  rost  insects,  micro-organisms  or  plants  of  the 
grass  family. ^ 


i  N., 

.  \  ’ 
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"reservoir"  which  consists  of  the  various  inorganic  sources  (the  soil, 
the  atmosphere,  the  lithosphere,  artificial  fertilizers,  etc.)  and 
sinks  (including  ultimate  consumers,  e.g.  man  or  livestock). 

Let  us  now  examine  the  matrix  elements  and  see  how  many  can  be 
determined  quickly.  The  various  elements  arc  as  follows: 


Matrix  Element 
A,|(i-^l) 

d.|  2  0  2) 

3  0-^3) 

a2i(2*-»l) 

^23(2-^3) 
^24  (2  “^4) 

<*32(3-*2) 

<l^l(4-»l) 

«'42(4-*2) 

<i43(4-*3) 

^44(4->4) 


Interaction 

I 

Reproduction,  growth,  decay  of  vegetation 
Artthropod  feeding  on  vegetation 

Vertebrate  feeding  on  vegetation 

Decay  of  dead  plant-feeders  and  metabollj:  wastes 

Reproduction  &  growth  of  plant  feeders 

Predaceous  or  parasitic  feeding  on  arthropods,  etc. 

Vertebrate  feeding  on  arthropods,  etc. 

Decay  of  dead  predators/parasites  and  metabolic  wastes 


fPredSt ion/paras  it  ism  on  other  predators/parasites 
^Reproduction  &  growth  of  predators/parasites 

Vertebrate  feeding  on  predators/paras iteb 

I 

Decay  of  dead  vertebrates  and  metabolic  was.:os 

I 


{Vertebrate  feeding  on  other  vertebrates 
[Reproduction  5  growth  of  vertebrates 


It  is  immediately  obvious  that  ■  ^32  ■  >  ^43!”  since 

the  corresponding  interactions  do  not  exist  in  the  particpiar  ecosystem 
being  considered.  Furthermore,  it  is  quite  reasonable  to  assune  that 
^21  ^  ^31  ^  *ln^^  'ew  arthropods  die  a  natural  dea  ;h;  most  are 
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cons'jmed  alive  sooner  or  later  by  other  arthropods  or  vertebrates.  To 
a  lesser  extent  this  Is  also  true  of  vertebrates,  but  among  the  domi¬ 
nant  groups  (birds  and  rodents)  death  by  starvation  Is  sufficiently 
frequent  to  warrant  attention.  We  thus  have 


r<«||  4,2  0 

0  d22  ^23 


'41 


0  433  434 

0  0  444 


Initial  state 


As  a  consequence  of  the  perturbation  discussed  earlier  (mammals  and 
birds  destroyed),  the  last  row  and  colunn  would  be  missing  In  the 
perturbed  matrix. 


a,,  a;,  •  \ 

0  «22  <*23  I  final  state 

K  °  233/ 


We  allow  for  the  possibility  of  a  non-tero  value  for  4j,.  Since  both 
matrices  are  time-dependent  arrfJ  tiie  system  is  open,  we  have  ten  Inde¬ 
pendent  4's  and  ten  phase-constants  to  determine  In  order  to  specify 
the  Initial-state  matrix  completely.  Six  9^s  and  six  phase-constants 
are  sufficient  to  specify  the  final-state  matrix. 


Possible  "Conservation  Laws" 

I.  All  phase  constants  remain  unchanged,  e.g. 

Tij  ■  Tij  15. 

assuming  meteorological  conditions  remain  the  same. 


i 


t 

] 

I 
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2.  On«  would  tend  to  expect  non-zero  diagonal  elements  to 
remain  unchanged,  since  the  rate  at  which  organisms  in  a  given  class 
reproduce,  grow  and/or  feed  on  each  other  Is  unlikely  to  depend  strongly 
on  what  happens  elsewhere  In  the  ecosystem.  Thus  (In  the  example), 

«*ll— a<22  **^2*  ^33”  *33 

3.  Conservation  of  total  ^  In  a  closed  system.  In  this  case 
equations  (l)  would  be  applicable  and 


17. 


In  an  open  system,  such  as  our  particular  example,  modified  approximate 
versions  of  (2)  and  (17)  may  still  be  available.  Examination  of  each 
individual  case  Is  necessary  In  the  chcsen  example  it  seems  reasonable 
to  set 


max  2^  4 1  m  max  Aj 
j-l  j-l 


18. 


which  expresses  the  likelihood  that  the  maximum  amount  of  protoplasm 
produced  In  the  Initial  and  final  systems  would  be  very  nearly  the 


same. 


4.  functional  relationships  for  ^-transfers  can  be  used  If 
they  are  sufficiently  simple  and  well  established.  In  the  c:iosen 
example.  It  seems  reasonable  to  assume  that  the  rate  of  arthropod 
feeding  on  plants  is  essentially  proportional  only  to  arthropod  popula¬ 
tion,  with  the  same  proportionality  before  and  after  the  perturbation, 
'thus,  .  " 

^12^*2^  *12^  '9- 

which  incidentally  impi  lesT|2  »  A2  20. 

Admittedly  the  insect  population  depends  indirectly  on  the  amount 
of  food  available,  but  this  may  be  thought  of  as  a  second-order  rela- 
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ttonshtp  and  is  seldom  Iimitln9.*  Initial  conditions  (number  of  over¬ 
wintering  larvae,  etc.)  and  weather  factors  are,  of  course,  averaged 
out  in  the  model . 

We  also  conjecture  that  the  rate  of  arthropod  predation  and 
parasitism  on  other  arthropods  is  approximately  proportional  to  preda- 
tori(>a  rasite  populat  ion,  and  that  the  ratios  are  conserved: 

^23/^3  “  ^23^^3 


This  in  turn  implies  ■^23  -^^3  ”, 

A  variety  of  equations  describing  functional  interrelationships 
may  be  suggested  for  more  complex  ecosystems,  but  it  must  be  recognized 
that  the  basis  for  many  of  them  is  likely  to  be  theoretical,  if  not 

hypothetical,  to  a  degree.  To  the  extent  that  the  number  of  possible 

•a 

equations  involving  the  <^'s  exceeds  the  number  needed  to  determine 
completely  the  perturbed  matrix,  the  more  uncertain  relationships  can 
be  tested  against  one  another.  This  is  not  our  present  purpose, 
however,  since  the  chosen  example  is  too  limited  in  scope  to  make  such  a 
sophisticated  procedure  seem  likely  to  yield  worthwhile  results.  It 
must  be  remembered,  on  the  other  hand,  that  in  a  more  complex  model 
there  will  be  correspondingly  more  unknowns  and  the  question  if  use¬ 
fully  interpreting  some  degree  of  over  l  iterminat ion  is  moot. 

For  our  example  we  can  use  the  thi-be  equations  (I^),  plus  (id), 

(19)  and  (21),  to  fix  the  six  independent  S's.  Thus  the  perturbed 
-matrix  Is -  - - -  -  --  -  -  ^  -  - 


*0d>jti  et  al.'*  point  out  that  Insects  seldom  approach  maximal 
feeding,  whereas  rodents  often  do.  Presumably  this  reflects  stricter 
supervision  of  insect  populations  by  predators,  and  great  sensitivity 
to  weather  conditions. 
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Mil  ‘*12  o\ 

j  0  ^22  ^23 1 

\sii\  0  e(33/ 

From  the  earlier  discussion  of  equation  (9)  It  will  be  recalled  that 
there  Is  a  linear  homogeneous  set  of  equations  involving  the  time- 
independent  amplitudes  3  defined  by 

-^j)]  -S'jk 

and  the  amplitudes  Aj.  The  conditional  for  solubility  is  Cet  3f  «  0, 
which  becomes 

0  cos  -r,j))  -i]  pjzCiAcos  6>„(T22  -T23))  - 

A  cos  ^qCTii  ~  ^3))  *  0  ^*^12  ^23  ^31  Acos  <u>o(  'fi\  -  A3) 

23. 

In  analogy  to  equation  (8)  we  have 

0  -[a,,  iAcostO^(T,,  -  ■r,2)  -  i]%  ♦  ^31  i  Aeos  t^^(r3,  .  ^3)23 

0  -  iS|2  A,  -*.[322  iAcos  <Oq(  J22  -  T23)  -  ^  % 

•*  25. 

0  -  ia23%  *^33  iAcos  (■ar33  -  X3)  -  j]^3 


Substitute  (19)  In  (25)  and  (21)  In  (26),  yielding 
^Acos  ^0(^22  *  ^23)] 

*2  ■  ’  ^33  iAcos  dL>o(r33  -  A3)]] 

Whence,  using  (19)  again,  we  obtain 

^'2  *  ^  f  ■  ^33  «Acos  iJo(T33  -  ^  3^ 
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Equation  (I8)  calls  for  maximizing  the  sum  of  several  sinusoidal 
functions  with  different  phases.  This  would  involve  a  good  deal  of 
rather  tedious  labor,  and  for  the  sake  of  simplicity  we  replace  (id) 
by  a  relation  involving  amplitudes  only: 

K  ■  Ai  +  A2  +  A3  +  Ai,  »  A|  +  A2  +  A3 
We  have  introduced  the  constant  K  for  .'“nvenlence.  From  (27),  (28) 
and  (30)  we  obtain 


A3  -  K  -  4 


-  ^22  i-^cos  ^o(V22  ‘  ^23^] 


^1  -  5133  i/i.cos  Oo(  ar33  -  ^  3^^ 


Aj 

^23  L.’ 

Now  insert  (28)  and  (31)  in  (25)  to  get  ^23' 


30. 


31. 


Finally,  using  (23),  (29)  and  (32),  we  can  determine  <5*3],  We  shall 
not  bother  to  exhibit  ^23  ^31  explicitly.  This  completes  the 

formal  analysis  of  the  model  ecosystem 

Perhaps  the  most  interesting  result  is  equation  (29),  which  deter* 
mines  the  rate  of  insect  feeding  on  vegetation  in  the  perturbed  case. 
Qualitatively,  the  conclusions  can  be  stated:  arthropod  feeding  in  the 
perturbed  case  will  be  maximal  for  systems  where  arthropod  predation 
and  parasitism  in  the  unperturbed  case  were  minimal,  as  indicated  by 
small  populations  (A3)  and  small  growth  potential  (^33)-  The  greater 
the  difference  between  the  phases  ^^3  and  ^3,  the  greater  the  effect 
of  Insect  feeding  in  the  perturbed  case. 


|8 


It  would  be  Interesting  and  Instructive  to  carry  through  a  graphi¬ 
cal  sna4ysis  of  the  equations.  It  vould  also  be  desirable  to  use  the 
iterative  step-function  approach  on  the  same  model  ecosystem— perhaps 
assyming  four  discrete  Increments  of  three  months  each.*  However,  the 
mathematical  complexity  of  the  Job  is  rather  higher  than  was  the  case 
with  our  chosen  example,  and  it  might  pay  to  program  the  iteration  for 
a  computer  and  go  to  a  much  larger  number  of  iterations.  Here  the  ob¬ 
jective  would  probably  have  to  be  to  follow  the  actual  time  behavior 
of  a  system  where  a  sufficient  number  of  matrix  elements  was  known 
explicitly  as  functions  of  time.  Predictions  could  be  made  of  the 
behavior  of  the  remainder.  Our  more  ambitious  objective  of  predicting 
the  response  of  an  ecosystem  to  perturbation  seems  beyond  the  scope 
of  the  discrete  model  (simply  because  less  is  assumed  and  more  has  to 
be  calculated) . 

Sources  of  Data 

To  complete  the  illustration  we  shall  outline  s  number  of  possible 
data  sources  or  experiments.  Care  must  be  exercised  to  choose  suffi¬ 
ciently  many  Independent  sources,  e.g.  it  would  be  Inconsistent  to 
measure  two  quantities  related  by  a  "conservation  equation"  separately, 
while  using  the  equation. 

It  Is  Important  to  keep  In  mind  that  Is  n^t  simply  the  5  (e.g.  nitrogen) 
content  of  the  standing  crop  of  the  k^^  class.  No  allowance  need  be 
made  in  the  continuous  model  for  that  portion  of  the  standing  crop 

^Actually  there  is  no  particular  advantage  In  choosing  equal  time 
periods.  The  technique  permits  great  flexibility  In  the  choice  of 
Increments. 
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which  is  "consumed"  by  the  ecosystem,  but  the  portion  excreted  or  used 
up  in  metabolic  processes  must  be  counted  as  though  it  were  still  phy¬ 
sically  present.  On  the  other  band.  In  the  discrete  model  the  potent  la  I 
crop  during  a  given  increment  of  time  must  also  be  considered,  and  if  the 
time  increments  are  considerable,  the  necessary  corrections  to  allow  for 
turnover  may  be  sizable. 

In  the  case  of  grain  crops  the  potential  crop  during  a  growing 
season  is  probably  of  the  order  of  20-25%  larger  than  the  actual  crop,** 
but  there  Is  no  N-e.xcret ion.*  In  the  case  of  short  1  ife-cyc I e  organisms 
the  difference  may  be  much  greater.  Odum  et  aj..®  estimate  the  biomass  of 
the  potential  crop  of  grasshoppers  in  an  old-field  ecosystem  to  be  of  the 
order  of  10  times  the  observev'  maximum  standing  crop  ot  peak  periods  (an 
estimate  which  is,  admittedly,  uncertain  at  best).  Unfortunately,  census 
counts  are  only  suitable  for  measuring  the  standing  crop,  whence,  from 
the  point  of  view  of  collecting  useful  data,  the  continuous  model  must 
be  preferred,  Cven  here  the  difficulties  are  severe.  One  vexatious 
problem  is  to  deduce  the  total  population  from  the  size  of  a  series  of 
samples,  since  an  absolute  count  is  generally  infeasible.  If  not  impos¬ 
sible,  for  arthropods  and  other  small  organisms.  Some  recent  work  Indi¬ 
cates  that  real  progress  is  being  made  in  overcoming  these  difficulties,^ 
but  results  are  prel iminary  so  far.  Horeover  metabolic  losses  of  nitro¬ 
gen  In  excreta  are  difficult  to  measure  in  roanv  cases  (e.g.  arthropods) 
and  this  fact  introduces  further  uncertainties. 

To  Make  the  problem  concrete,  iet  us  specify  a  particular  eco¬ 
system  for  vdileh  existing  data  are  reasonably  concrete,  namely  a  homo¬ 
geneous  grain  fatmlng  area  such  as  Champaign  county,  Illinois.  Much 
*lf  we  were  tracing  the  H2O  economy,  transpiration  would  be  Impor- 

* 
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basic  information  is  available  in  publications  of  the  College  of  Agri¬ 
culture  of  the  University  of  Illinois.^  For  the  particular  county 
in  question  the  average  crop  mix  Is  43.2%  corn,  31.2%  soybeans,  and 
much  smaller  percentages  in  oats,  wheat,  forage  and  woodland.  The 
"growing  season"-- measured  arbitrarily  from  the  mean  date  of  the  last 
spring  freeze  to  the  mean  date  of  the  first  fall  freeze— embraces  the 
peri'od  April  25  to  October  20  180  days).  ^ 

The  Influence  of  bioci  imatological  factors  on  phenology  has  been 

0 

Investigated  comprehensively  by  the  American  Institute  of  Crop  Ecology.” 

It  has  been  found  that  growth  and  maturity  can  be  indexed  quite  accu¬ 
rately  in  terms  of  the  number  of  day-degrees*  elapsed  between  different 
stages.  Extrapolating  some  of  these  results  to  the  soil  and  climate  of 
Champaign  county  (using  the  local  conditions  to  "normal  I ze")  should 
make  It  possible  to  estimate  with  reasonable  accuracy,  apart  from 
the  contribution  arising  frem  decay  organisms.  Unless  it  is  desired 
to  focus  attention  on  the  decay  process  explicitly,  almost  any  reason¬ 
able  assumption  can  be  made  regarding  the  latter,  providing  only  that 
consistency  is  maintained.  In  fact,  in  an  open  model  the  decay  process 
can  be  lumped  with  the  external  interactions  (the  reservoir)  vnd  i 

ignored  altogether,  providing  the  interpretation  of  ^||  is  consistent. 

Other  data  of  a  rather  cruder  sort  are  also  available  In  the 
literature.  For  example,  there  exist  estimates  of  the  net  primary 

production  (in  grams)  of  various  crops  under  different  conditions.^  . 

For  corn,  in  particular,  Illinois  approaches  the  ecological  optimum 

'  I 

with  yields  of  1011  grams  per  square  meter  per  year,  or  5.6  grams  per 

I 

square  meter  per  day  during  the  l80-day  growing  season. These  '  ^  1 

figures  must,  of  course,  be  translated  into  N-content  but  many  analyses  I 

■  I  —  II  ■  _  IW  .!■  .11  . .  I  ^  II  .  II  . .  ,  I 

^Measured  above  a  kO  F.  base.  ^ 
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of  the  chemical  composition  of  various  plant  parts,  for  various  crops, 
ar>^  avallabte.^^  These  integrals  average  out  the  t  ime*dependence  of  the 
quantities  in  question  but  are  useful  as  a  means  of  double  checking. 

By  working  backwards  (If  the  sinusoidal  time  behavior  is  assumed)  the 
amplitude  Ai  can  be  deduced  independent  I y~-at  least  for  corn.  Still 
uther  types  of  data  are  available,  such  as  detailed  measurements  of 
N-uptake  from  soils  by  various  crops.  Here,  of  course,  alloivance 
must  be  made  for  contributions  by  N-fixing  bacteria  living  in  the  soil.^^ 

If  a  forest  >yere  being  considered  the  estimates  would  be  .nore 
difficult  to  make.  A  recent  technique  has  been  discussed  wherein  CO2 
production- is  measured  and  used  as  an  Index  of  metabolism.*^  Estimates 
of  the  standing  crop  In  term*  of  canopy  foliage,  understory  foliage, 
woody  parts,  etc.  have  been  made. 

Standing  crops  of  arthropods  and  other  animals  are  much  more 

difficult  to  estimate  and  subject  to  much  greater  fluctuations.  For 

farmlands  the  best  source  of  primary  data  on  grazing  insects  is  probably 

lb 

the  Cooperative  Economic  Insect  Report.  Beneficial  insects  are  also 
occasionally  reported,  but  evidently  less  thoroughly. 

Obtaining  the  N-content  of  body  material  plus  wastes  i'  also  a 
problem,  especially  for  arthropods,  since  protein  percentage  alone  is 
not  a  reliable  guide  in  this  case.*  Unfortunately  the  typical  arthro¬ 
pod  exoskeleton  (integunent,  or  chitin)  has  a  molecular  structure  in 
which  certain  N-containing  amino  acids  are  linked  tightly  to  long- 
chain  molecules  similar  to  cellulose.  Both  the  N-content  and,  more 
important,  the  absolute  quantity  of  Integiment  vary  widely  from  species 


*As  a  rule  of  thirib,  on  a  dry-weight  basis  protein  averages  I6X 
nitrogen. 
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to  ^neclei,  making  extrapolations  from  the  few  well  known  cases’^  diffi¬ 
cult.  In  spite  of  the  problems,  however.  It  should  be  possible  to  esti¬ 
mate  A2  with  reasonable  confidence,  though  there  are  probably  Insufficient 
census  data  to  obtain  A^. 

Losses  of  crops  by  Insect  feeding  are  also  estimated  from  time  to 
time  by  the  USOA  in  terms  of  a  percentage  of  the  monetary  value  of  the 
national  harvest. T->  the  extent  that  national  averages  are  appli¬ 
cable  to  the  chosen  geographical  area,  and  to  the  extent  that  damage 
to  the  seed  or  fru’t  portion  of  the  crop  reflects  damage  to  other  por¬ 
tions,  some' sort  of  estimate  for  0()2  might  be  derived.  Since  the  seed 
or  fruit  constitutes  a  small  percent  of  the  total  plant  protoplasm 
(averaging  7%  for  herbaceous  annuals  and  grasses) — which,  however, 
varies  considerably  from  species  to  species—and  since  arthropods  as  a 
class  do  not  feed  on  all  plant  parts  equally  (though  their  exact  degree 
of  preference  for  the  more  succulent  parts  would  be  difficult  to  esti¬ 
mate  quantitatively),  there  are  very  many  uncertainties  in  this  type 
of  calculation.  Laboratory  or  greenhouse  experiments  may  be  possible 
which  would  aid  In  these  determinations,  especially  as  regards  the 
composition  of  body  wastes,  although  It  can  be  extraordlnar ■ ly  d I f f Icul t 
to  simulate  field  conditions.  Another  possible  avenue  is  the  use  of 
radio- isotopes  as  tracers.  A  number  of  projects  using  tracer  techniques 
are  under  way  at  present,  at  Oak  Ridge  National  Laboratory  and  at  the 
University  of  Minnesota,  among  others.  It  seems  possible  that  some  of 
their  data  may  be  applicable  to  the  present  problem,  even  though  the 
experiments  were  carried  out  with  different  objectives  In  view.  Tracer 
techniques  may  offer  the  only  practical  possibility  of  determining  Aj, 
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and  with  a  good  deal  of  refinement  it  is  conceivable  that  an  approxima¬ 
tion  for  0(23  could  be  obtained. 

The  two  dominant  vertebrate  populations  are  likely  to  be  the 
birds  and  rodents.  Bird  populations  are  strongly  influenced  by  the 
presence  or  absence  of  cover  (woodlands  or  large  trees)  and  in  Cham- 

\ 

paign  county,  where  only  1.3%  of  the  land  is  wooded,  suitable  nesting 

sites  are  somewhat  limited.  An  average  of  3  birds  per  acre  has  been 

'  / 

estimated  for  Illinois  farmlands,  inc I ud I ng  horned  larks,  starlings, 
crows,  nuthatches,  bluejays,  cardinals  and  a  number  of  r  in,;r  species. 

Census  counts  are  made  by  numerous  local  birdwatching  clubs  and  kept 
in  the  archives  of  the  National  Audubon  Society.  The  Bureau  of  Biolo¬ 
gical  Survey  is  another,  perha.'S  more  rel  iable,  source,  both  for  birds 

18 

and  other  animals.  Counts  of  rodents,  shrews,  voles,  moles,  bats, 
skunks,  etc.  are  probably  an  order  of  magnitude  less  complete  and  may 
be  altogether  inadequate.  However,  direct  experimental  verification 
in  geographically  limited  areas  should  not  be  unduly  difficult,  since 
the  ground-living  animals  (unlike  birds)  may  be  considered  essentially 
immobile.  N-oontent  of  warn  -blooded  animals  is  not  difficult  to  com¬ 
pute  from  the  weight,  but  metabolic  losses  (excreta)  are  dom>r:nt. 

An  average  U  oz.  meadowlark  is  likely  to  consune  about  2U  times  its 
own  weight  in  food  each  year,  and  a  f  oz.  shrew  may  consume  two  or  three 
times  its  own  weight  per  day.  Detailed  correlations  between  body 

weight  and  metabolic  rate  for  warm-blooded  animals  are  available  In  1 

5 

the  literature*^,  whence  extrapolations  from  one  species  or  one  stage 

in  the  life  cycle  to  another  should  be  within  the  scope  of  possibility.  | 

The  diets  of  at  least  some  common  species  have  been  studied  with  consid-  i 

\ 

I 
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erable  care.  Very  detailed  analyses  of  the  seasonal  variation  of  the 
stomach  contents  of  the  comnon  crow  and  starling  have  been  published, 
and  others  doubtless  exist.  By  thorough  analysis  of  available  d^ta, 
and  perhaps  a  minimal  addition  of  new  data.  It  should  be  possible  to 
estimate  A^,  d^\,  rt|4.  A24>  <(34  444. 

In  Summary,  It  appears  feasible  to  obtain  more  or  less  direct 
experimental  evidence  which  would  be  capable  of  yielding  numerical 
values  for  most  of  the  A's  and  fli's.  However,  in  view  of  the  uncer¬ 
tainties  involved,  one  would  feel  relatively  little  confidence  in  deter¬ 
minations  of  ^i2t  <^23’  ^22’  ^33>  ^2  ^3*  functional  rela¬ 

tionships  relating  some  of  the  elements,  e.g. 

Ai  +  A2  •  A3  +  A4  ■  K  30. 

^3  “  •^"*3 

analogous  to  (19)  and  (20)  must  be  used,  or  uncertain  numerical  estimates 
must  be  accepted.  In  practice,  the  best  procedure  would  probably  be 

to  Iterate  the  solutions  by  repeating  the  calculations  using  different  » 

sets  of  assumptions  until  the  maximum  degree  of  self-consistency  Is 
achieved.  Ultimateiy,  of  course,  better  techniques  of  expert:' f*nta I 
verification  must  be  found. 
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